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1 Introduction 



In recent times the interest in understanding hard processes involving nuclei, like electron nucleus 
collisions (e + A) and hadron nucleus collisions {h + A), has rapidly increased. The relativistic 
heavy ion collider (RHIC) will provide a new tool for experimental studies of h + A and A + A 
physics. This experimental advance provides a challenge to theorists for developing a description 
of such reactions strictly in terms of the underlying fundamental field theory, i.e. quantum 
chromodynamics (QCD). 

The usual way to proceed in QCD is to use factorization theorems ^ that permit a rigorous 
treatment of hard processes. In this approach, incalculable soft contributions are factorized into 
universal, process independent distribution functions. Measurement of the distribution functions 
in one process then allows the prediction of cross sections for other processes. An example is the 
Drell-Yan (DY) process where two hadrons scatter to produce a lepton pair of large invariant 
mass Q^. If is large enough the factorization theorems state, that the inclusive Drell Yan cross 
section can be expressed as the product of a calculable hard part and universal twist-2 parton 
distribution functions that can be measured, for instance in deeply inelastic lepton nucleon 
scattering. Such universal distribution functions have been successfully extracted from existing 
data. The natural question arises how this picture is changed for hadron-nucleus collisions. 
Differences between the twist-2 distribution of a single nucleon and those of a nucleus of mass 
number A > 1 were first discovered by the European Muon Collaboration. This phenomenon is 
known as the EMC effect at large values of the Bjorken scaling variable x and as shadowing at 
lower values of x 0. 

However, strong interaction dynamics cannot solely be described in terms of twist-2 structure 
functions. These functions are simple one-particle correlators and have the interpretation of a 
probability to find a parton in the hadron with certain momentum fraction x. These distribution 
functions are measured when the hard probe scatters off a single parton in the nucleon. Of course, 
the probe may scatter off more than one parton in the target. Such processes are related to multi- 
parton correlators and describe corrections to the single-parton scattering process. In general 
they are suppressed by powers of the large scale, and are called contributions of higher twist. 

If we consider processes involving large nuclei it is quite obvious that multiple scattering will 
be enhanced and even dominate in the limit A ^ oo. The classical example for enhancement 
of multiple scattering in h + A processes is the Cronin effect, i.e. the anomalous A-dependence 
of the Drell-Yan cross section at large values of the transverse momentum of the lepton pair . 
Instead of a scaling of the cross section with the volume ~ A as expected in the twist-2 case, 
one observes an A-dependence parametrized by a^"^"^ = A"a^~^^ with a 4/3 This is a clear 
signal for multiple scattering and shows that higher-twist multiple scattering is grossly enhanced 
at large momentum transfer. 

The literature contains a wide variety of approaches dealing with multiple scattering, see 
e.g. 0, but a rigorous formulation in the language of perturbative QCD was attempted only 
recently. In a series of papers, Luo, Qiu and Sterman (LQS) and recently also Guo, showed how 
to extend the factorization theorem beyond leading twist to get a handle on double scattering in 
terms of perturbative QCD |^. Their twist-4 matrix elements are correlators of two partons 
in the nucleus and can describe scattering off partons from different nucleons. This approach 
predicts a stronger dependence on the size of the nucleus than the linear scaling with A of single 
scattering. LQS argue that — due to the geometrical structure of the matrix elements — some 
of these twist-4 nuclear matrix elements pick up an additional factor A^^^, compensating for 
some of the inherent suppression of higher twist corrections. This nuclear enhancement is able 
to describe the observed scaling with A^^'^ in certain kinematical regions. 
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Figure 1: Single- (a) and double-scattering (b) schematically for a hadron h colliding with a 
nucleus A producing a fi'^fi~ pair. One additional unobserved parton is radiated in the leading 
contribution to high transverse momentum q±. 



In this paper we discuss nuclear enhanced double scattering contributions to DY pair pro- 
duction off nuclear targets 

h{P2) + A{P,A)^l+l-+X (1) 

in the framework of the LQS approach. It was already shown some time ago by Guo p|, ^ that 
higher-twist effects in the process (|l]) are enhanced at large transverse momentum. Extending 
the calculations of Guo 0, we have calculated the full kinematical dependence of the cross 
section da/dQ^ dq'j_ dydfl including the angular distribution of the lepton pair and evaluated it, 
in particular, for p + A scattering at RHIC. Our most important results were already published 
as a letter In this contribution we present the complete analytical results together with 

a brief review of the LQS approach and some technical details. We also give some numerical 
results not contained in and an extended discussion. We will argue that nuclear enhanced 
double scattering corrections cannot be neglected at RHIC. 

In Fig. Ill we show typical Feynman diagrams contributing to the DY process. Single scattering 
means that, e.g., a quark from the single hadron and an antiquark from the nucleus annihilate 
into a virtual photon in a hard process. To obtain large transverse momentum an additional 
unobserved parton, here a gluon, has always to be radiated. Double scattering in our example 
implies that the quark can pick up an additional gluon from the nucleus. The generalized twist- 
4 factorization theorem of LQS states that two distinct scattering reactions occur for double 
scattering. One is referred to as double-hard (DH) scattering. Here the parton from the single 
hadron undergoes two subsequent independent hard scattering reactions in the target. The 
other one is called soft-hard (SH) scattering. Here the parton from the single hadron interacts 
with the collective soft gluon field in the nucleus (it can e.g. be thought of as feeling the color- 
magnetic Lorentz force) before producing the virtual photon via the final hard scattering. At 
large transverse momentum interference of these two mechanisms can be neglected. 

The double-hard reaction resembles the classical double scattering picture, while soft-hard 
scattering may have an analogy in plasma physics, where a hard scattering event can be ac- 
companied by initial- or final-state interaction with the mean field. The first process might be 
accounted for by event generators in simulations of heavy-ion reactions whereas the implemen- 
tation of soft-hard scattering has yet to be formulated in a fully consistent manner (see [|1^ for 



a conceptual discussion). 

It would be helpful to find observables with vanishing or at least very small contributions 
from single scattering. The measurement of such observables would permit the identification of 
twist-4 effects at larger values of where perturbation theory is more reliable and thus provide 
a clearer evidence for the occurrence of double scattering. We find these observables in the 
angular distribution of DY pairs, namely the helicity amplitude Wa and the Lam- Tung relation 
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Waa = [|T2|) both introduced in section 2. The Lam- Tung relation is analogous to the 

Callan-Gross relation in deep- inelastic electron nucleon scattering (DIS). Furthermore it would 
be helpful to find observables where double-hard and soft-hard contributions can be disentangled. 
These can also be found in the angular distribution of DY pairs. While SH scattering leads to 
non-trivial angular patterns, DH scattering does not. Any anomalous A-dependent angular 
distribution therefore would be a clear sign of non trivial twist-4 correlations. 

Although we will focus on h + A collisions, our work has been primarily motivated by the 
fact that higher-twist effects or multiple scattering should be most prominent in nucleus-nucleus 
{A + A) collisions at high energy, as they will be routinely studied at RHIC. Two microscopic, 
Q CD- inspired descriptions of nucleus-nucleus collisions at RHIC energies and beyond have been 
proposed in recent years. In one, the parton cascade picture [O, O, [13], the collision is modeled 



as a sequence of incoherent binary scattering events among partons ("minijets"), which can be 
cast in a probabilistic formulation as a multiple scattering process. In the other one, the random 
mean-field model |jl6[ the collision is described as a nonlinear interaction among coherent color 
fields carried by the colliding nuclei. We hope that the study of the interplay between soft-hard 
and multiple-hard processes in a simpler framework, such as h + A interactions, will help to 
elucidate the competition between these processes in A -|- A collisions. 

The determination of multiparton correlation functions in nuclei is also an important step 
towards the application of QCD to nuclear physics in general. The generic enhancement of 
higher-twist effects in reactions involving nuclei facilitates their extraction and give the chance 
for a deeper theoretical understanding. They encode the essential aspects of the difference in the 
quark-gluon structure between isolated nucleons and nuclei. From the practical point of view a 
better understanding of such medium effects is necessary to be able to interpret the so called 
"hard probes" of high-energy heavy-ion collisions in a convincing manner. 

The unintegrated Drell-Yan cross section is actually the classic example for a typical two-scale 
process. Both and q±, the mass and transverse momentum of the photon, are detected. The 
perturbative treatment then gives rise to logs of the type \og^ {Q'^ / q'^) which can become large if 
the scales differ substantially. In such a case they have to be resummed. Resumming large logs 
in the presence of additional soft gluons contained in the twist-4 correlators is a delicate task 
that could not be solved yet. We therefore follow a less ambitious procedure and require that 
both scales are of the same order, thus reducing the original two-scale problem to an effective 
one-scale problem. In addition the LQS formalism in its present form requires qj_ to be not 
much smaller than in order to neglect interference terms that are not calculated up to now 
and which would introduce additional twist-4 matrix elements. Recently Guo, Qiu and Zhang 
proposed a method to describe the region of low transverse momentum and how to match it to 
the perturbative calculation at high q^ . 



Finally we would like to comment on the experimental situation. There exist data for angular 
coefficients onn + A scattering, e.g. from the NAlO |T^ and E615 experiments which are not 
yet completely understood in terms of conventional Q CD-treatment. These experiments measure 
pair masses Q ^ 4 . . . 8.5 GeV and transverse momenta q± < 3 GeV which is certainly at the 
limit of where we could still trust our calculation. We will focus on future RHIC perspectives 
where we could hopefully reach higher transverse momenta. 

The paper is organized as follows: In section 2 we give the general definition of the Drell-Yan 
kinematics, we introduce our observables — helicity amplitudes and angular coefficients — and 
explain possible choices of the photon rest frame. In section 3 we remind the reader of the 
known twist-2 results for the helicity amplitudes. Section 4 is devoted to the calculation of the 
double-scattering contribution. We explain some essential features of the factorization procedure 
of LQS. In particular we discuss the distinction between double-hard and soft-hard scattering 
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and give an introduction of how to calculate them practically. Finally section 5 contains our 
numerical results along with a discussion and we close with a summary of our final conclusions. 
In the appendices we have collected some technical details referred to in the text. 



2 General Definitions 
2.1 Kinematics 

Let a nucleus with momentum per nucleon collide with a hadron with momentum P^. We 
study lepton pair decay of virtual photons with momentum and mass = q^. In the hadron 
center-of-mass (cm.) frame we take 

P{' = i(v^,0,0,V5), P2'^ = ^(v^,0,0,-v^) (2) 

where ^/S is the total center-of-mass energy. The differential cross section for Drell-Yan pair 
production is then given by 

The solid angle Q is parametrized by the polar and azimuthal angles 6 and of one decay lepton 
in the rest frame of the virtual photon. We have the freedom to choose the axes in the photon 
rest frame conveniently. We will discuss this point after eq. ([T6|). 

The hadronic tensor follows the standard definition and is given by 

W,, = j d^xe^^-^' (P1P2I j^(x)j,(0) IP1P2) . (4) 

The leptonic tensor is 

U^^ = 2{r,r,+^ll)-g^^^Q\ (5) 

where li and I2 are the momenta of the leptons. 

We rewrite the g-dependence of the cross section in terms of rapidity transverse momentum 
and mass Q of the virtual photon in the hadron cm. frame. Introducing the hadronic 

Mandelstam invariants 

S = (P1 + P2)' 

T = (Pi-g)2 (6) 
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+ 1 fQ^-U 



the rapidity y of the photon is 

Likewise we can express the transverse momentum of the lepton pair in the hadron cm. system 
by 



and we can rewrite the differential cross section as: 

do" a^^vr 



em 



dQ^dqldydn ASQ^ ^'^ (27r 



L,.W^''7K-U (9) 
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2.2 Helicity amplitudes and their frame dependence 

To discuss the angular distribution of the lepton pair it is convenient to introduce hehcity am- 
phtudes Way defined as 

Way = e^^h^:'>W^- (10) 

where e!r\ a E {±1,0} is a set of polarization vectors of the virtual photon with respect to the 
axes defined by the spatial unit vectors X, Y, Z in its rest frame: 

e^^^ = ^{TX-zY)„ ef = Z, (11) 



Following |£0l it is advantageous to introduce a complete set of structure functions 

= Wi^i 
= Wo,o 

Wa = ^ (W^i,o + Wo,i) ^^^^ 

which are referred to as the transverse, longitudinal, spin fiip, and double spin fiip helicity 
amplitudes. Due to parity conservation other combinations are related to these by Way = 
W-a-a'- In terms of these helicity structure functions the cross section can be written 

as 

do- al^ 



dQ2 dql dy dQ 647r35Q2 

(wtl (l + cos^^) + Wl Q - ^cos^^^ + iyA(sin2^cos0) + iyAA (sin^ ^ cos 20)^ . (13) 

The angular dependence can be integrated out which reduces the cross section to 

dcr aL / 16 \ ,^ ^, 

—71 Wtl (14) 



em 



dQ2 dql dy 647r3Sg2 

with Wtl = Wt + |Wl = ~^Wl^. The ratio of the angular dependent and angular integrated 
cross sections is usually parametrized in two standard forms [0, |21| : 

I6n ( dQ2 dg2 Ay dQ \ ,^ n a 3 o r\ * ■ c^n i ^2 . v ^ ^ , 

da — 1 = (1 + cos^ 9) + Aq[-- - cos^ 9] + Ai sm 29 cos (/> + ^ sm'' 9 cos 2<^ 

dQ2 dg^ dy 



A + 3 



^1 + A cos^ 9 + fi sin 29 cos ^^^^ ^ ' 



(15) 



in terms of two sets of angular coefficients 



^0 = TTT- M = — — A2 = , (16) 

Wtl Wtl Wtl 

2-3Ao 2A 2A2 .... 

A = 11 = v = :-• 17 



6 




At this point it is necessary to comment on the proper choice of coordinate axes in the photon 
rest system with respect to which the polar and azimuthal angles 9 and are defined. Clearly 
the helicity amplitudes and angular coefficients depend on the chosen frame. Two prominent 
examples are the CoUins-Soper frame (CS) and the Gottfried- Jackson frame (GJ) which are 
discussed in detail in |^0|]. When we boost from the hadron cm. frame to a lepton cm. frame, 
the collinearity of the hadron momenta Pi and P2 is lost and they span a plane which we identify 
with the X-Z-plane. We are still free to fix Z within this plane. In the Collins-Soper frame we 
choose Z to bisect the angle between Pi and — P2, see Fig. 0. The angle between Z and Pi 
and between Z and — P2 is called 7cs In the Gottfried- Jackson frame Z lies in the direction of 
Pi. We have introduced the the angle 7gj between Pi and P2 for short notations in that case. 
Obviously we have 27cs + 7gj = ^r. See Appendix 7.1 for more details. Also note that for g_L = 
both frames are equal and is no longer defined. 

In order to obtain the results in both frames in a straightforward way we use the projectors 
given by Mirkes ||2^ 



TL 



L12 



flU 



-9 



-'1 



E1E2 



^h2 



p^pu 

El 



where Ei and E2 are the energies of the hadrons in the lepton cm. frame. In terms of hadronic 
Mandelstam variables and the photon momentum these are 

Q^-T ^ e-y 1 ^,^, , ^ Q''-U 

2Q 



El 



S{Q' + ql), E2 



S{Q' + ql). 



(19) 



a G 



2Q 2g V 2Q 2Q 

Contracting the hadronic tensor with these projectors gives amplitudes 
{TL, LI, L2, L12} which are just linear combinations of the helicity amplitudes in the different 
frames. These can be found by linear transformations 



( 



V<2 



M, 



(/3) 



TlI 
Tl2 
V Tl12 / 



(20) 
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with (3 G {CS,GJ}. The transformation matrices M(^) follow from simple geometry and are 
given in Appendix 7.1 for the CS and GJ frames. 



3 Leading twist results 

To make this paper self contained we give the leading order results for the angular coefficients. 
Due to twist-2 factorization the hadronic tensor is given by a convolution of a perturbatively 
calculable partonic tensor and two twist-2 distribution functions: 

W,. = g^elY, ^H^,t\^u^2)fa/Ai^i)Mm^5is + t + u-Q'). (21) 

a,b 

The sum runs over partons a and b in the nucleus A and the single hadron h respectively and 
^1, ^2 are their momentum fractions. We denote the strong coupling by g and is the fractional 
electric charge of the quark flavour coupling to the virtual photon. The /(^) correspond to 
twist-2 parton distribution functions deflned in the usual way as 



MO = p^j ^P^e^^^-y- \ (P| mi^Vq{y-) IP) 
fg/hiO 
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1 27r 



(22) 



for quarks and gluons, respectively in a hadron moving with momentum P+ along the light cone 
|ll]. q and F denote quark flelds and gluon fleld strength tensors and V is the path-ordered 
exponential 

P = Pexp jzc/^ dx"A+(a;-)| (23) 

which makes the deflnition gauge invariant. In the following we will not write this path-ordered 
exponential explicitly, but it is assumed to be inserted between all factors in a product of 
operators referring to different space-time points. In the = gauge this term will disappear 
but we will always use covariant gauge. Note that we also suppressed the dependence of the 
distribution function on the factorization scale /i^. 

The partonic tensor H^^'^ represents the hard part of the process. We have introduced 
partonic Mandelstam variables by 

t = (^,P,-qf = ^,(T-Q') + Q' (24) 
u = {^2P2-q? = UU-Q') + Q'. 



The (5-function in eq. (^) stems from the on-shell condition for the unobserved emitted parton 
with momentum /. It can be used for trivial integration of ^2 via 

We can decompose the partonic tensor in a way completely analogous to the decomposition of 
the hadronic tensor W^i, in terms of helicity amplitudes. Therefore we deflne the four projections 
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Figure 3: Twist-2 contributions to order as', annihilation (a), and Compton processes (b and c). 
The blobs L and R stand for all possible tree-level diagrams. 



i^TL^; -f^L^^5 H^'', which are obtained in the same way as in eqs. ( p^Ojp^ . This allows us to 



write hadronic tensor and helicity amplitudes as 



W^ = 8 



1 



B 



6 UQ'-T)-Q' 



fa/A{il)fb/h{i2)Hl^^{il, i; 



(26) 



with B 



-u 



6 



Hgi(T-Q^) 



g^l^T^ and a G {/iz/, TL, L, A, AA}. There are 



three kinds of partonic processes contributing to the leading-twist process at large transverse 
momentum: quark-antiquark annihilation (g + q) and Compton processes with a gluon from the 
single hadron (g -|- ^f) or the nucleus [g -\- q). They are shown in Fig. ^. 

The cross section for the leading twist DY pair production at large transverse momentum is 
For completeness we give the results in the CoUins-Soper frame in Appendix 

T2[. It can be 
or 2i/ — (1 — A) =0. This relation holds for 



well-known ^ 
7.2. 



ii^L is precisely the so-called Lam- Tung relation 



The relation -ffAA 
expressed in terms of angular coefficients as Aq 
all partonic subprocesses and is not specific for the Collins-Soper frame but holds in any photon 
rest frame as long as the Z-axis lies in the reaction plane. Corrections arise at next-to-leading 
order (NLO) H. Note that Ef^ picks up an additional sign compared with if^"*"" in contrast 
to the other helicity amplitudes. This is due to the inversion of the orientation of the angle 7cs 
if we interchange the roles of P\ and P2- For symmetric collisions, e. g. J9 + p, this leads to the 
leading- twist result Wt\ = 0. 



4 Twist-4 Factorisation 

4.1 Soft-hard and double- hard scattering 

For the calculation of the double scattering contributions we use the approach of Luo, Qiu and 
Sterman which is presented e.g. in We consider it worthwile to repeat some aspects here. 

In Fig. 1^ we illustrate the forward diagrams that need to be calculated for double scattering. 
The box denotes all possible partonic subprocesses. We have to introduce a cut to obtain the 
corresponding hadronic tensors. This is not only possible in the middle but we also have to 
include interference graphs which appear if the forward diagram is cut in an asymmetrical way. 
In Fig. ^ we give concrete examples for graphs which involve a quark-gluon matrix element in 
the nucleus and an antiquark coming from the single hadron. The first diagram is cut in the 
middle (M), the others are cut left (L) and right (R). 
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Figure 4: Collection of the schematic forward processes we want to consider for double scattering. 
Factorization allows separate treatment of twist-4 and twist-2 matrix elements (T and /) and 
hard part (H). 



The starting point for the calculation of the hadronic tensor is the forward scattering ampli- 
tude in fourth order of the strong coupling constant g: 

= ^ y d4xd^Xid"x2dSdx4e^^^ {PiP2\T[Jp{xs)Ux,)j^{0)j,{x) ^27) 

Mx2)Ja{x,)A''{xs)A-{x,)AP{x2)A'^{x^) + ...] \P,P2) 

where j (J) are electromagnetic (color) quark currents and A denotes gluon fields. T is the time- 
ordering operator. Dots indicate further contributions involving more gluon fields with gluon 
self interaction and less quark currents such that the power of g remains four. By application of 
Wick's theorem we obtain all diagrams which potentially contribute to the forward process. 




Figure 5: Examples for double-scattering diagrams with the partonic process qg + q cut in the 
middle (a), right (b) and left (c). Circles in Fig. (a) indicate poles from propagators which lead 
to finite gluon momentum (open) or soft gluon momentum (filled). 
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We classify diagrams by their partonic subprocess ab + c where a and b are partons from the 
nucleus and c comes from the single hadron. Our notation for the parton momenta is shown in 
Fig. ^. We concentrate first on qg + q which can be realized e.g. by the diagrams of Fig. ^. As 
a first step towards factorization we perform Fierz transformations in color and Dirac space for 
the fermion lines entering and leaving the blobs: 

QaiOHi^) = -\ [-^Jf^'^^^ + K'ti) +■■■■ (28) 

Here the are the SU (Nc) Gell-Mann matrices, and Nc is the number of colors. Latin characters 
are color indices, greek one are Dirac indices. We keep only Dirac structures which give leading 
contributions to the unpolarized Drell-Yan process. After Fourier transformations of the fields 
we end up with 

T-+^-- /^^^^^ /d^d^z.dS.dS. 

(27r)^5(^)(ri + r2 + K - I - q)e'^^^^e-'''ye'^''e-'^'''^ (Pi| q{0)j\{zi)Af{z3)A''{z^) \P,) 

1 1 

( - ^) {P2\ mYqiy) \P2) e?^? V/v^ 5Z sid)Z,,,- (29) 

We used the fact that forward matrix elements only depend on the relative coordinates. So, 
after shifting positions in each matrix element separately, we introduced a new set of space-time 
coordinates [y^zi^z^^z^). I is the momentum of the unobserved radiated parton. The index d 
runs over all diagrams of the class qg + q and S{d)^^^^,^^ contains all perturbative propagators 
and vertices of the hard part of diagram d, with 7^ (7^) inserted for the upper (lower) blob, and 
traces over color and Dirac indices. For the diagram of Fig. ^(a) we have e. g. 



cRS = 1 Tr 



'^^''>1 71? tFi \2 I • (~^7/3^ )-( 1 , ■ (-^7a^ 

{Ti + K — K' — qY + le (r2 + K'Y + le 



,R. ^(- h- ^) I ^(/i- k) 



(r2 + Ky + le (ri — qY + le 



P + ie 



(30) 



R and 5* are color indices of gluon fields. Note that we have coupled quark and gluon field 
operators in the nuclear matrix elements separately to color singlets. We have omitted a term 
proportional to {d^^^ + if^^^) {Pi \ q'y^t^qA^''A^" \Pi) which would induce color forces between 
quark and gluon in the nucleus. This would not allow to factorize quarks and gluons in color 
space to test the size of the nucleus. Such contributions are not expected to show nuclear 
enhancement, i.e. scaling by A^^^. 

Next we perform the light cone expansion, i.e. we use that the partons move almost coUinear 
to the hadron in the infinite momentum frame. We introduce the parton momentum fractions 
by setting ri = ^iPi, r2 = ^2-^2, K = xPi + Kj_ and K' = x'Pi + K'j_. r[ = ri + K — K' is fixed by 
momentum conservation in the upper forward matrix element. Note that we have to allow one 
parton (called b here) from the nucleus to be soft with essentially zero longitudinal momentum. 
Therefore we must take into account transverse momenta K± and K'j_ to get a contribution 
in higher orders of the coUinear expansion as we will see later. Now one can integrate the 
trivial degrees of freedom in eq. (^). To complete the factorization we perform a Sudakov 
decomposition of the 7-matrices in the forward matrix elements on the light cone and retain 
only the leading terms 7'' = -^7 , 7"^ = ^7~- Then we can introduce twist-2 distribution 
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functions fc/h{^2), given in eq. (22), for the lower blob. Finally we cut the forward diagram. 
This sets the radiated unobserved parton with momentum / on the mass shell and in the hard 
part all expressions right to the cut receive a complex conjugation. 

Hence after these first steps we have factorized the hadronic tensor in a straightforward way. 
We get 



J a « 



(31) 



f„,(ei, X, x', K^, K'^)H:'+'{^,, 6, X, x', K^, K'^). 



We have derived eq. ^ for our examples of Fig. ^ but this equation holds for all processes shown 
in Fig.|. We have separated the soft part involving the twist-2 distribution function fc/h and 
the nuclear matrix element Tab from the hard part H'^^'^. 

For our example qg + q we are left with a quark-gluon matrix element in the nucleus 



f X X' K< K',)- I ^^3 d'^^i ir+z7 JK+z- -iK , -iK'+z7 JK', 

i,g[i;i,x,x - I 2^ ^2^^3 ^2^^3 e e e e e 



(32) 



For other classes of processes we also need the quark-quark and gluon-gluon correlators Tgg and 
Tqq which are given by eq. (l3^ ) with the second line replaced by 

-1^ (Pil F-+iO)F^iz^)APiz^,zi)A^iz,,zi) \P,) for f,, 

^1^1 ^ (33) 

- {p,\miM^i)qi^3, 4)rQ{^4, 4) l^i) for Tqq. 

Again here the field operators of each single parton form color singlets separately. We have 
arranged the first two field operators referring to the nuclear parton a in the way they appear in 
the corresponding twist-2 structure function while the remaining two field operators of parton b 
are still in their original form. 

The hard part for the qg + g-process is 

HX%^i,^2,x,x',K^,K'^) = i^^^J-^(-^"/^)rK^E^(^)-w4„t- 

Hard parts for other processes are obtained analogously. The delta function in eq. (|3TD can be 
rewritten for the cases in which the diagrams are cut in the middle (M), right (R) or left (L) as 



(35) 



(M): 


m 


= (65 + T- 




-'6{^i + x- 6) 


(R): 




= (65 + T- 




-'6{^i+x-x' -Q 


(L): 




= (65 + T- 




- ^c) 



where 

^b = ^^s + T-Q^ = 61^^=0. (36) 

As pointed out in various papers 0, §] there are in principle two contributions to double 
scattering: There can be two hard scattering reactions of the parton c from the single hadron 
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or the outgoing unobserved parton on one side and partons a and b from the nucleus on the 
other side. This is called double-hard scattering. Furthermore there can occur processes where 
the parton c from the single hadron or the outgoing parton picks up one soft nuclear parton b 
in addition to one hard scattering off the nuclear parton a. This is called soft-hard scattering. 
These are not only qualitative notions but they can be rigorously derived from the equations 
given above. One can understand them in terms of poles in the hard part of the scattering 
tensor. In eq. (^Tj) we have four momentum integrations d^i d^2 dx' and four propagators. 
We will perform the ^i-integration by making use of the delta function The propagators 

provide poles with respect to x and x'. In general we expect two poles in x and two poles in 
x'. Application of the residue theorem will fix the momenta and it is easy to check that two 
poles give large momenta x = Xhard = x', respectively and the other ones give small momenta 
X = Xsoh = x' with Xhard > a;soft -C 1 as long as q_L ~ Q. 

To be more precise the pole structure we expect for our example in Fig. |^(a) is 

W., ~ f dx- / dx'- ^ -, (37) 

J (a; - Xsoft + ^e)(a; - Xhard + «e) i (a;' - Xsoft - «e)(a;' - a;hard - ^e) 

where L and R stand for the numerators of the amplitudes left and right of the cut depending 
on X and x'. The residue theorem implies 



L(x — Xsoft) L(^X — Xhard) 



•^soft -^hard -^soft -^hard 



R{x' = Xsoft) R{x' = Xhard) 



•^soft -^hard -^soft -^hard 



(3^ 



Soft-hard and double-hard scattering correspond to the quadratic terms proportional to L(x = 
Xsoft)R{x' = Xsoft) and L{x = Xhard)-R(a;' = Xhard) respectively. We follow the arguments of 
reference ^ and neglect the remaining mixed terms of eq. (|3^) which are interference terms 
between soft and hard rescattering. This is permitted if q± is not too small. In the case qj_ <^ 
however we obtain Xgoit ~ a^hard- The interference will then be important and eventually spoil 
nuclear enhancement due to the negative signs in eq. (|38| ) [§. Since the integrated Drell-Yan 
cross section does not exhibit nuclear enhancement, we even expect a suppression effect at low 
values of transverse momenta 0. 

Most diagrams we naively get from Wick's theorem for double scattering show a different 
behaviour. For them the propagators do not provide four but only three or two poles — but at 
least one in each variable x and x'. This implies that only soft-hard, only double-hard or neither 
of these contributions can occur. For example it is easy to understand from this analysis that 
only symmetrically (M) cut diagrams can contribute to double-hard scattering. Consider the 
diagram of Fig. ^(b). The propagator on the right-hand side has denominator [{r[ — q)'^—ie], from 
which X and x' dependences cancel out. On the left-hand side, the propagator at the bottom 
of the diagram contributes a denominator (r2 + + ie = {x + ^ + ie). This is a soft 

pole which fixes x to Xsoft = ^ 1- The remaining two propagators on the left-hand side 

fix x' to either a soft or a hard value. Therefore we have two contributions: (i) x = x' is small. 
This means that parton b is soft — in Fig. ^(b) this is the gluon from the nucleus. Together 
with the production of the virtual photon which is always a hard process, this constitutes a soft- 
hard double scattering process, (ii) x is soft and x' is hard. This is what we called a soft-hard 
interference because parton b has different momenta on the left and the right side of the diagram, 
i.e. momentum is transferred between partons a and b. We omit this contribution as explained 
above. We conclude that there is no double-hard scattering contribution from the diagram in 
Fig. |^(b) since x cannot be hard. On the other hand the diagram in Fig. Ka) is an example in 
which all four poles occur, indicated by small circles. 
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We now discuss another technical point which has some relevance for nuclear enhancement. 
Consider again eq. ( p?]) and the example of Fig. |^(a) and note that the numerators L and R 
contain exponential functions e*^^i ^3 g*(5i)-^)^'i ^1 q-^^'Pi . To apply the residue theorem 
we have to close the integration contours. For — zjf > we have to close the x-integration 
in the upper complex x-plane with no poles encircled and therefore vanishing integral. On the 
other hand for 2;;^ — < we must close in the lower complex x-plane with soft and hard pole 
encircled. Therefore these poles always introduce additional theta functions in coordinate space, 
in our example 0(Z]~ — 2;^) from the x integration and 0(— zj) from the x' integral. The presence 
of these G functions in the matrix elements is an essential ingredient of the proof of the following 
rule about the cancellation of final state interactions for soft-hard scattering: If at least one soft 
parton from the twist-4 matrix element is directly coupled to the outgoing unobserved parton 
line these graphs do not contribute to soft-hard scattering. These graphs exhibit more than one 
possibility to set the cut and the nuclear enhancement cancels between the contributions with 
different cuts due to the particular arrangement of 9-functions in that case [0]. 

Hence, after a careful discussion of the pole structure of each diagram the number of dia- 
grams we are finally forced to evaluate can be reduced. Figs. ^ and ^ show the diagrams that 
contribute to double-hard and soft-hard scattering, respectively. We marked the poles used to 
fix momenta. Also we omitted soft-hard contributions with soft quarks since they should be 
suppressed compared to their counterparts with soft gluons. 




a) -^M b) -^M c) -^M d) -^M 



Figure 6: Graphs contributing to double-hard scattering: qg + q^ qg + g, qq + g and qq + q. There 
are additional processes with structures similar to the last one, e.g. qq + q. This is discussed in 
the next section. The blobs (L) and (R) stand for all possible tree-level diagramms. Propagators 
providing poles which are used to fix hard momenta are indicated by small open circles. 




a) -^M b) -^M c) -^M 



Figure 7: Graphs contributing to soft-hard scattering: qg + q, qg + g and gg + q. The blobs 
(L) and (R) stand for all possible tree-level diagramms. Only symmetrically cut diagrams (M) 
are shown. Corresponding right (R) and left (L) cut diagrams have to be added. Propagators 
providing poles which are used to fix soft momenta are indicated by small filled circles. 



14 



4.2 Calculation of the double-hard processes 

First we investigate the case of double-hard scattering. Since all parton momenta are hard we 
can work in the coUinear limit and set = = K'^ immediately. We perform the integrations 
over X and x' only taking into account contributions from hard poles. These poles are marked 
with circles in Fig. |[ They fix momenta to the values x = Xh = x' and = Xa with 

Xa = (39) 

Note that Xa is determined by kinematical variables only whereas Xh depends on the integration 
variable ^2- Two additional 0-functions in the light cone space-time coordinates show up as 
explained above. 

The field operators of the second parton in the nuclear matrix element are selected in such 
a way that we get the usual twist-2 light cone operator q'j^q for quarks and field strength 
tensors F^^F^^ for gluons. In the latter case we have to turn two factors of into light cone 
derivatives {ixh)~^d^- , {—ixh)~^d^'- and perform partial integrations to let them act on the gluon 
fields in the matrix element. Hard gluons carry only physical polarizations, this allows us to set 
d"^AP{z^)d^A''{z4) = F'^-^{z^)F+^{z^)l{-gP'' + nPn" + nPfi") with the usual light cone basis 
vectors nP = Pf/P^, nP = P^/P2- 

We are now able to give the final result for double-hard scattering. For a partonic subprocess 
of the type a6 + c we have 

= ^'"^^1'!"/'^ /J 1^ 6^ + T - g2 /c/Hfe)g(x.,a:.)^i7r"^+-(a:.,a;.,6) (41) 

with a G {/ii/, TL, L, A, AA}. Note that parton a is always a quark or antiquark for double- 
hard scattering and Cg denotes its electric charge. As in the helicity amplitudes which arise from 
contractions of the hadronic tensor W^^y defined in eq. |1^, here we use the same contractions 
for the hard partonic part H^^, defining amplitudes i^TL, -f^L, -f^A and -f^AA depending on the 
frame. The hard parts for the relevant diagrams can be written as 

H^^^'^''-^ = c^-JA{h- Mp«(/i- ^)7m hiu] 6=.. la'-^r 
^2"'"'^' = [(/i- /i)R<^P /^p«(/i- M a=.. \ a'-^r (42) 

^/ ^ X=x'=Xf^ ^ 

^5 ^ X=x'=Xfi ^ 

Lpa and R^is indicate the sum over all possible tree level diagrams for the blobs at the left and 
the right side of the corresponding diagram in Fig. p. Each diagram provides an individual 
color factor which is indicated symbolically by the factor C. For partonic processes with three 
participating fermions we have to distinguish between different possible combinations of quark 
and antiquarks and different flavours. In principle there are contributions written in our short 
notation as qq + q, qq + q and qq + q. Each of these contributes four diagrams contained in 
Fig. H(d). In addition, there are processes of the type qq + q', qq' + q and qq' + q', which are each 
represented by a single diagram. In our notation g is a quark or antiquark, q the corresponding 
antiparticle and q' denotes a quark or antiquark of a different flavour than q. As an example we 
give the hard part for the qq + q subprocess: 
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- Tr [(/i- ^)7, A^7a(/i- ^)7/. /i7.] (n + - g)-'(r2 + K)'^ + 

+ (K^l)} ^,=^, (43) 

Note the relative sign that is due to the different number of fermion loops in the corresponding 
diagrams. 

The T^i^{xa,Xh) are the universal nuclear matrix elements for double-hard scattering intro- 
duced in P, They depend on the momentum fractions Xa, Xh of both hard partons from the 
nucleus and are given by 

rS''(^a,x,) = ^^ I dz4-^^e(^r-z3')e(-^r)e^^"^^^^i"e^^'^^^^(^3--r) 

(Pil Mm^h^Qih) \Pi) 

(Pil g(0)7+g(2;r)g(2;3-)7+g(^4 ) \Pi) 
(Pil g(0)7+g(^r)g(C)7+g(^3") ■ 

We are now ready to explain the origin of nuclear enhancement in these matrix elements. As 
mentioned above the fields of both partons are coupled separately to colour singlets. This 
allows to probe distances up to the diameter of a large nucleus. Each spatial integration of the 
matrix element arises from a Fourier transformation and we expect it to be accompanied by a 
rapidly oscillating exponential function (since — > oo). At first glance, this seems to rule 
out coherence over the whole nucleus. However only and — in eq. (0) are restricted 
by phase factors while z^ + z^ is not. Hence, we can set / d{z^ + z^) to be proportional to 
the nuclear radius Ra ~ tqA^^^ and obtain an additional factor of A^/^. Here ro ~ 1.1 fm is the 
radius of a nucleon. 

To find a reasonable model for the double-hard matrix elements imagine the insertion of a 
complete set of states Ylip \P) {P\ I {'^P^V) = 1 where V is the volume of the nucleus. Neglecting 
for a moment the 0-function and the integration over z^ + z^ we can compare these definitions 
with the common definitions of twist-2 structure functions in eq. (p2|). Assuming that the sum 
over states is saturated by the lowest lying nucleon state double hard matrix elements reduce to 
a product of two twist-2 structure functions times a dimensional factor which scales like A^^^ 

) = CA^^^fa/A{Xa)fb/A{Xh). (45) 

fa/A is the nuclear parton distribution per nucleon in a nucleus with mass number A and proton 
number Z , N = A — Z , given by 

Z N 

fa/A{x) = -jfa/p{x) + —fa/n{x). (46) 

In a primitive model one can simply set the remaining free integration with respect to z^ + z^ 
equal to Ra and obtains (7 = 3/ (SvrrQ) ^ 0.005 GeV^. Although the factorization of the double- 
hard matrix element into a product of two twist-2 distribution functions and the scaling by A'^^^ 
seems to provide a reasonable model, the choice of the normalization factor C is a hot matter of 
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discussion, see section 5. The 0-functions provide an ordering of both single scattering events 
on the hght cone. They simply express the fact that scattering off parton b takes place before 
the scattering off parton a. This corresponds to < on the left side and to < on the 
right side of each diagram. 

Now we give the final results for the hard parts. For the evaluation it is important also 
to include ghosts, whenever one can close gluon lines through the non-perturbative blobs. For 
H^^ = — iif^^^'^, which enters the angular integrated cross section, we reproduce the results 
of [§]. We have listed them in Appendix 7.3 for all subprocesses. The helicity amplitude Wtl is 
independent of the frame. Due to the special form of all diagrams contributing to double-hard 
scattering we easily see that the contractions P^^H^^ and P^^2H]i^ in eq- (11) always vanish, 
and for the last contraction we get 



L2 t.. {^Q2 _Tf(^Q2 _uf TL (g2 + ^2 )2 TL V I 

for all processes. For the Collins-Soper frame this leads to the simple relations 



W^^= sin^csW^i?H =jy^W^^ (48) 

sin7cscos7csW^TT = -q^W^tl" (49) 

W2l= il^r- ^ (50) 

Immediately one observes that the Lam- Tung relation holds in the CS frame and all other 
frames considered here. For the Gottfried- Jackson frame the results get particularly simple. The 
double-hard contributions for deviations from the simple 1 -|- cos^ 9 behaviour vanish completely: 

W^'' = {), Wf^ = Q, W2l = {). (51) 

At first glance this seems to be rather surprising but it is a direct consequence of the special 
form of all double-hard contributions and therefore has a simple explanation. 

The double hard process resembles the classical double scattering picture. To understand 
this we consider again the process depicted in Fig. |I]. The antiquark from the hadron undergoes 
a first scattering with a hard gluon from the nucleus. By subsequently radiating a gluon with 
large transverse momentum the antiquark returns to the mass shell, but now with transverse 
momentum q±. The second scattering then corresponds to the lowest-order annihilation process 
of two massless quarks. Therefore the double-hard process can be factored as the product of two 
single scattering cross sections, supported by our model that the double-hard matrix element can 
in good approximation be written as a product of two twist 2 distribution functions. This process 
naturally fulfills the Lam- Tung relation and depends only on the structure function WpL ~ ^ll- 
Other helicity amplitudes differ only by kinematical factors, i.e. can be obtained from WpL in 
the Gottfried- Jackson frame by simple rotations. In this probabilistic picture scaling by A"^/^ 
appears to be very natural. The cross section of the first scattering obviously scales by the 
volume ~ A of the nucleus. Afterwards the projectile parton travels further through the nucleus 
and has to penetrate the nuclear matter. The probability whether a second interaction takes 
place on the way to the boundary of the nucleus scales then by the nuclear radius ~ A^^^. 

4.3 Calculation of the soft-hard processes 

For soft-hard scattering we have to keep K_[^ in the expansion of the gluon momenta. The limit 
= would lead to if = which gives no contribution to physical double scattering but just 
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a contribution of the eikonal phase. Note that we are always using covariant gauge where these 
factors are present. Starting from equation (^) we do a colhnear expansion 

of the terms depending on intrinsic transverse momentum 



We have used K\ = K\ which allows us to write 



^ = K^—^ + K',^ ^. (54) 

^ dirt ^dK^ ^ dK'^^ 



The mixed terms of second order proportional to K^K'^^ give the leading contribution to soft 
hard scattering. The factor K^K'j_'^ can be converted into partial derivatives with respect to , 
. A partial integration and an expansion = A'^Pf /P^ + . . . for soft gluon fields in the 
nuclear matrix elements turns them into field strength tensors. The K± and K'j_ integrals can 
then be carried out. Note that having soft quarks with transverse momentum and converting the 
corresponding Kj_ into covariant derivatives acting on the quark fields would lead to operators 
contributing only at the level of twist-6 0. 

We perform the pole integrations for the symmetrically cut graphs (M) in Figs.0 and the 
corresponding asymmetrically cut ones (R) and (L). The residues fix the momentum fractions of 
the soft gluons to x = Xg = x' while the hard partons carry ,^1 = Xf, for the symmetric diagrams 
and = Xc for the asymmetric diagrams and 

= % (55) 



Xh — C,b Xs — t C I /n2 t C ^ 



We have introduced the notation fc^ = —K^ > 0. When we take a forward diagram for soft-hard 
scattering then the sum of differently cut diagrams (M), (R) and (L) can be divided into two 
terms. This procedure is described in detail in |^. In the first term, phase factors and B- 
functions arrange in such a way that this term is bounded in space-time by geometrical arguments 
and cannot pick up an additional factor A^^^. This contribution shows no nuclear enhancement. 
The second term is not bounded and only depends on the hard part of the symmetric diagram 
(M). We therefore neglect the first term and arrive at 

' d6 



.. r™(.,)^ff™*'(.„.„6) (58) 



Xk 



4 ^ dK} dKf 



1 



Kj^=0 Xb 



where a G {/iz/, TL, L, A, AA}. The hard part for the qg + q process shown in Fig. 0(a) is 



(59) 
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and Ri^/s stand for the tree-level diagrams on the left- and the right-hand side of the diagram 
in Fig. |^(a) indicated by the blobs. C symbolizes the color factor of each single diagram. Hard 
parts for qg + g and gg + q subprocesses are obtained in an analogous fashion. 

The universal matrix elements appearing for soft-hard scattering depend only on the mo- 
mentum fraction Xb of the hard parton and are given by 

1 r r\ - r\ - ^^^^ 

for scattering off a quark-gluon and a gluon-gluon pair in the nucleus, respectively. We can 
repeat the arguments given below eq. for the nuclear enhancement. Since the longitudinal 
momentum fraction of the gluons is essentially zero, we even encounter only one phase factor 
gixi,p^ z-^ contrary to the case of double hard scattering discussed before. One therefore might 
expect that soft-hard matrix elements lead to an enhancement factor proportional to the square 
of the nuclear radius A^^^ since only one out of three integrations is limited by phase factors. 
However, the fields at positions and are correlated by color forces. The anticipated absence 
of long-range color correlations in a nucleus implies that integration over z^ — z^ is still restricted 
to a region of the size of a nucleon. 

To build a model for the soft hard matrix elements one again can compare eq. (^) with the 
definition of the usual twist-2 matrix elements eq. (^2]). Neglecting for a moment the G-functions 
we see that the twist-4 pieces differ only by an additional Xb independent renormalization. As- 
suming that this can be factored into a term proportional to the nuclear radius and a (possibly 
scale dependent) constant 

~ / dz:,^F-^z^)F^ M^) (61) 

the Xb dependence of the twist-4 matrix element can be modeled by the usual twist-2 structure 
function. 

Tj(x) = X'Alfa/A{x) (62) 

Eq. ( |6lD can be interpreted as the collective color-magnetic Lorentz force experienced by the 
projectile parton while traveling through nuclear matter p. 

The only i^'^-dependence of the hard part will be implicit in Xb and Xs- Therefore we can 
rewrite the derivatives in the form 

' d6 1 



^ (x„42jj ^^^S + T~Q^y^ 



(63) 



+ T!^{xc)-{H'J''^'^i{x^,^: 
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where H{xc,^2) = H{xc,Xs = 0,^2) and 6) = ^^H{xc,Xs = 0,^2)- Note that in g] 

the last term containing a derivative with respect to Xg is absent since it turns out that 
does not depend exphcitly on Xg. We give the results for the hard parts in the Collins-Soper 
frame in Appendix 7.4. There are some subtleties in removing the unphysical degrees of freedom 
for processes of the class qg + g. There we have soft quarks from the nucleus with dominating 
component ^4+ and for the hard gluons from the single hadron we have to take only the transverse 
polarizations. For Wf? reproduce the results of ^j. Unlike for double-hard scattering the 
amplitudes Wl^, and cannot be obtained in a trivial way from PVtl- Related to that 
is an obvious violation of the Lam- Tung relation. 



5 Numerical Results 

In order to cover different energies S we calculated helicity amplitudes and angular coefficients 
for two different settings: for 252 GeV tt" on tungsten as in the Fermilab E615 experiment 



T9|| and for future p + A collider experiments at RHIC energies with 250 GeV protons and 
A X 100 GeV large nuclei. The main task before doing the numerics is to insert a suitable model 
for the twist-4 matrix elements in the nucleus. The models introduced in the last sections for 
the double-hard matrix elements 

T^''{Xa,Xt,) = CAlf,/A{Xa)fb/A{Xh). (64) 

and for the soft-hard matrix elements 

T!^ix) = X'Alf,/Ax) (65) 

are already given by LQS 0, If. We tried to argue above that their dependence on the parton 
momentum fractions and their scaling behaviour with respect to A seem to be quite natural. 

Nevertheless the correct normalization of the models is not yet clear. The hope is that C 
and A are universal constants. In the early days of the model a value of C = 0.072 GeV^ was 
taken from theory p], ^ while was taken from comparison with experiment. This was done 
by analyzing dijet photoproduction on nuclei 0, ^ leading to a value A^ = 0.05 — 0.1 GeV^. 
Later an analysis of g_L-broadening in hadron-nucleon Drell-Yan processes implied a value of 
A2 = 0.01 GeVl Obviously there is a discrepancy, which may either point to a strong scale 
dependence or questions the universality of the soft-hard nuclear matrix elements. Note that the 
soft interaction besides the hard process in soft-hard double scattering is a final state interaction 
for dijet photoproduction but an initial state interaction for the Drell-Yan process. 



In their latest work Guo, Qiu and Zhang |17| used the fact that there is a formal connection 



between soft-hard and double-hard matrix elements via the limit 

hm XgT°"(x„ xg) = T^f{x,) (66) 

which leads with the models of eqs. (|H^),([H3|) to 



C ^ . (67) 

Xgg\Xg}\xg^Q 

They take Xgg{xg)\xgaiQ ~ 3 as the gluon momentum density at Xg — >■ 0. Taking into account 
the later value of A^ this leads to a normalization factor C which is an order of magnitude below 
the former theoretical value. But it is more consistent with the value of C ~ 0.005 GeV^ which 
we obtained in section 4 from a naive normalization argument. 
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Figure 8: The angular coefficients A, /i and v (from left to right) as functions of in the 
Gottfried- Jackson frame for 252 GeV tt" on tungsten at Q = 5 GeV, y = 0. The dashed line is 
the result for single-scattering and the thick line is given by the full calculation with single- and 
double-scattering. We also show results of the E615 experiment at low transverse momentum. 

In another, more geometrical, model one can estimate the two-parton area density that can 
be seen by the projectile parton in the nucleus by the product of two single-parton area densities. 
Integration of the two-parton area density over the effective transverse area of the nucleus yields 
exactly the right-hand side of eq. ( |6lD with a constant C = 9/(87rrQ) ^ 0.01 GeV^. Obviously 
C is presently basically unknown and a true improvement can only be expected if comparisons 
of a larger number of observables with more precise experimental data are available. 

For soft-hard scattering we use the smaller value of = 0.01 GeV^ obtained from Drell-Yan 
data to maintain consistency. For double-hard scattering we have chosen the maximal value 
C = 0.072 GeV^ to give a feeling how big the effects of double scattering might be. We take 
the twist-2 distribution functions from the CTEQ3M parameterization for the nucleons |^ and 
from the recent GRS fits for the pion ||2^. Where nuclear effects on twist-2 distribution functions 
are taken into account the EKS98 parametrization |29| is used. 



Let us ffist address the question how this approach works with the pion data from E615. 
In Fig. P we show our results for the angular coefficients A, fi and u in the Gottfried- Jackson 
frame for the process 7r~ + W at E615 energy. We also give experimental data from the E615 
collaboration. We can fairly well describe the data for A at low transverse momentum but fail 
for fi and u. It seems to be an essential feature that corrections due to nuclear enhanced twist-4 
matrix elements can be quite large for the helicity amplitudes but the effects largely cancel for the 
angular coefficients which are ratios thereof. Therefore changes of the angular coefficients in the 
range q± ^ Q are small compared with the pure twist-2 calculation. Experimental values have 




-3-2-1012 -3-2-1012 -3 -2 -1012 



Figure 9: The helicity amplitudes wtl, "U^l and for p + Au collisions aX S = 10^ GeV, 
(5 = 5 GeV and g_L = 4 GeV. We show the twist-2 contribution (solid lines), double-hard 
(dashed lines), and soft-hard (dashed-dotted lines). Soft-hard is scaled up by a factor of 10. We 
give results with (thin lines) and without (thick lines) EKS98 modifications for nuclear parton 
distributions. 
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Figure 10: Violation of the Lam- Tung relation. Difference of angular coefficients Aq and A2 is 
shown. Calculation are with (solid line) and without (dashed line) EKS98 corrections to nuclear 
parton distributions. 



large error bars above q± = 2 GeV. We already pointed out that data at q±/Q ^ 1 would provide 
a much better test for the twist-4 calculation because of the absence of soft-hard interference 
terms and the absence of double logarithmic corrections of the type log^(Q^/g^) which have 
to be resummed. We have to conclude that nuclear enhanced double scattering contributions 
are not able to yield substantial improvement for the old pion data. Deviations of the n — W 
data from the naive perturbative Q CD-predictions might be explained by other approaches like 
the Berger-Brodsky mechanism which takes into account higher twist contributions of the pion 
wave function ||30|, in terms of non-trivial spin correlations in the QCD- vacuum |31|] or via 



contributions of the chiral-odd T-odd distribution function hf |32[ . 

We turn to p + Au collisions at RHIC with S = 10^ GeV^ per nucleon pair. We give helicity 
amplitudes instead of angular coefficients since these provide a clearer measure of nuclear en- 
hanced higher twist effects as explained above. To make numerical results of helicity amplitudes 
Wa, a € {TL, L, A, AA} more meaningful, we introduce the prefactor given in eq. (0), 

to convert the quantities into differential cross sections. In Fig. |] we give single scattering, 
double-hard and soft-hard results for the helicity amplitudes wtl, "M^tl and wa as functions of 



rapidity y in the Collins-Soper frame. One can read off the main results already given in |T0 
First the nuclear enhanced twist-4 contribution is absolutely important at RHIC energies. Even 
at quite large values of Q = 5 GeV and q± = 4: GeV used for the calculation its effect can exceed 
the leading-twist result for the parameters used. Second the double-hard contribution is by far 




0.01 

0.009 
0.008 
0.007 
0.006 
0.006 
0.004 
0.003 
0.002 
0.001 
0.0 



P+.4u 


b) 











Figure 11: Contributions of different subprocesses to wtl for double-hard (a) and soft-hard (b) 
scattering. Double-hard subprocesses are qg + q (solid line), qg + g (dashed line), qq + g (dotted 
line) and qq + q plus all pure fermionic processes (dash-dotted line). Soft-hard processes are 
qg + q (solid line), qg + g (dashed line) and gg + q (dotted line). 
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Figure 12: wtl as a function of transverse momentum q±_ at fixed Q = 'b GeV, y = Q (left) and 
as a function of photon mass Q at = 4 GeV and y = (right). Curves are shown for single 
scattering (solid line), double-hard (dashed line) and soft-hard scattering (dotted line). 



larger than the soft-hard part in the kinematical range under consideration. This can be seen for 
all helicity amplitudes in Fig. |. Even if we use smaller values of C proposed above double-hard 
scattering would reach the same order of magnitude as single scattering. 

We also give results with EKS98 modifications to nuclear parton distributions. These para- 
metrizations take into account shadowing, antishadowing and other nuclear effects on parton 
distributions, first described by the European Muon Collaboration. The numerical effect of 
these modifications on the single scattering results are small, but they turn out to be quite large 
for double-hard scattering. Note that it is not clear whether these nuclear effects for twist-2 
parton distributions should be carried over to the models for twist-4 matrix elements. We give 
results for modified and unmodified distributions to indicate the further uncertainty of the twist- 
4 matrix elements. The sensitivity of double-hard scattering to these EKS98 modifications is 
simply related to the fact that they involve a product of two nuclear parton distributions. The 
general question about the correct small-a; behaviour of the twist-4 nuclear matrix elements is 
very interesting since it has a big influence at RHIC energies. 

We return to the question of the Lam- Tung relation. We already saw analytically that double- 
hard scattering respects this relation while soft-hard does not. In Fig. |T0| we show the difference 
Aq — A2. It only picks up contributions from soft-hard scattering and therefore this quantity is 
very small. That means that violation of the Lam- Tung sum rule is a tiny effect for the process 
under consideration. 

The characteristic feature of the rapidity distribution is the peak of the double-hard distribu- 
tion at negative rapidities (i.e. for the photon emitted in direction of the initial single hadron). 
The reason for this is that negative rapidity prefers small Bjorken-x for partons from the nu- 
cleus. Parton distributions at small parton momentum fractions are large particularly for gluons 
and they enter quadratically for double-hard scattering. In Fig. ^ we separated the different 
partonic subprocesses contributing to double-hard and soft-hard scattering to enable a deeper 
analysis. One can see that double-hard processes with a quark and a gluon from the nucleus 
dominate over those with two quarks at negative rapidities (i.e. small Xh)- On the other hand 
qg + q is more important than qg + g in this region since momentum fraction ^2 of the parton 
from the single hadron is large and quarks dominate over gluons in that case. All the rapidity 
spectra can be understood in this way. 

Fig. 12 shows dependence of w^l on photon mass Q and transverse momentum g^. It can 
be seen that soft-hard and double-hard contributions as expected for higher-twist processes fall 
off more rapidly with Q and q± than (twist-2) single-scattering. 
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6 Summary 



We have calculated nuclear enhanced double-scattering contributions to Drell-Yan pair produc- 
tion in hadron-nucleus collisions using the framework of Luo, Qiu and Sterman. We found 
that the nuclear enhancement at RHIC energies and large transverse momentum makes double 
scattering a contribution that is equally important as single scattering. Measurements of an- 
gular coefficients and helicity amplitudes can distinguish double-hard and soft-hard scattering. 
Double-hard scattering has the interpretation of two sequential hard scatterings of onshell par- 
tons and is therefore related to probabilistic pictures of multiple scattering that are widely used 
in the phenomenology of heavy-ion reactions. In this spirit further comparision with alternative 
approaches to multiple scattering ^ would be very interesting. 



Double-hard scattering respects the Lam- Tung relation and, in contrast to soft-hard scatter- 
ing, does not show non-trivial deviations from the 1 + cos^ 9 angular distribution. Numerically 
double-hard scattering dominates over the soft-hard contribution. Soft-hard scattering violates 
the Lam- Tung relation but this violation is an effect of order 10^'^ for the model used. 

Double-hard scattering can be modeled as the product of two twist-2 structure functions 
which give the probability for finding two partons in the nucleus. The soft-hard matrix elements 
on the other hand encodes essential non-perturbative properties of the nucleus, describing the 
infiuence of the collective color field on the transverse momentum of the produced lepton pair. 
Inserting reasonable models for the corresponding matrix elements we predict, that collective 
color effects are negligible compared to the effects of subsequent independent hard scatterings in 
the kinematical range considered here. However the absolute normalization of the twist-4 matrix 
elements in these models is not yet clear. 

The double-hard scattering shows an interesting asymmetry in the rapidity distribution of the 
Drell-Yan pair. The pair is predominantly produced in the projectile, i.e. single hadron rapidity 
region. Although the absolute counting rates for the DY pairs at large transverse momentum 
are not high [|10| the predicted effect fits perfectly well in the PHENIX detector acceptance |]33 |. 



This asymmetry therefore should be observable. 

An analysis of Wl and W/\, particularly in the Gottfried- Jackson frame, as well as the 
measurement of M^tl (i-e. the angular integrated cross section) as a function of rapidity can check 
whether double-hard scattering is the dominant contribution. A possible violation of the Lam- 
Tung sum rule provides a further possibility to check the picture of soft-hard and double-hard 
scattering in general. If our predictions with respect to the ratio of soft-hard and double-hard 
scattering are confirmed by experiment it would provide support for the description of heavy 
ion reactions in terms of incoherent multiple scattering of free partons. Note that the prediction 
of a trivial frame dependence of double-hard scattering and the conservation of the Lam- Tung 
relation are independent of the models for the twist-4 distribution functions and provide good 
tests for the entire LQS approach. 
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7.1 Transformation of amplitudes to different frames 



Here we list the matrices M(^) defined in eq. (|20|) that map the invariant projections of the 
hadronic tensor to the set of hehcity amplitudes (Wtl, Wl, Wa, Waa)- They are 



cs 



GJ 



' 2 





1 

l\ 



- 



\2 





1 





1 



4 cos^ 7CS 



1 



4sin7cs cos7cs 
1+cos^ 7CS 
S cos2 7CS sin^ 7cs 





1 

cot 7GJ 

1+cos^ 7GJ 
2 sin^ 7GJ 









4 cos^ 7CS 

4sin7cs cos7cs 

1+cos^ 7CS 
8 cos2 7CS sin^ 7cs 






1 



sin 7GJ 



2sin7Gj 
COS7GJ 
sin^ 7GJ 







4cos2 7CS 



1-3 cos^ 7CS 
I cos2 7CS sin^ 7cs , 



(69) 



for the Collins-Soper and the Gottfried- Jackson frame, respectively. 7cs denotes the angle be- 
tween Pi and the Z axis in the photon rest frame and 7gj is the angle between Pi and P2. 
Obviously 27cs + 7gj = and in terms of hadronic Mandelstam variables we have 



COS7CS 



COS7GJ 



(Q2_T)(g2_f/) 
(Q2_T)(g2_f/)' 



sm7cs 



sm7Gj 



(Q2_T)(Q2_f/) 



(70) 



2Q2^ 



'(Q2_T)(Q2_f/) 



(g2_T)(g2_f/) 



1. 
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7.2 Single-scattering results 

Here we list the results for the hard parts Ha of single-scattering helicity amplitudes in the 
Collins-Soper frame. They are defined in eq. (|26|). We distinguish annihilation and Compton 
processes. 

H^^l^ = C,^{{Q'-tf + {Q'-uf) 



g2 _ t Q^-u 



tu \Q'^-t Q'^-u 
^TL --C2 - 



(g2-t)(g2_^) 

g^2(g2_t)2_(Q2_^)2 



E'^^ = -C 



fji+9 _ r 

^ St (g2-t)(g2-u) 



2 

2 „A2 



(g2-s)2 + (g^-M) 



su 

2 



^ _^ t (g^ + g)^ + (g^ - 
L (g2-t)(g2-M) 



g2t2(g2-M)2_(g2_t)2 



09+9 _ f 

A 2V (g2_i)(g2_^) 

Ci = Cp/Nc = 4:/9 and C2 = 1/{2N^ = 1/6 are essentially the color factors. 



(71) 
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7.3 Double-hard results 



In this section we give the hard parts i^xL of the projections for double-hard processes. 

They do not depend on the frame. 



27 V 65 Q^-TJ Q ^2S + T-Q^ 

TL or- \ ^1 ^ ' 



TTDH,qq+g _ ^ ( ' ~ j 

27 V 65 ^ 65 + T - g2 y ^ 6 (g2 _ 3^)^ 
^DH,..-+, ^ (65)^ + (65 + T - g^)^ (g^-Tf + (65 + T-g^ ^^' 

27 (g2-T)2 ^ (65)2 

4 (e25 + T-g2)2 



g2-T 65 + T-gv 16 (65)2 

65 + r-g2 X l(^^5)2^(^^^_^y_g2)2 



+ 



81 65(g2-T) 



rrDHM+-, ^ 1 [ (65)^ + (65 + T - g^)^ (g^ _ T)^ + 

27 (g2-T)2 (65 + T-g2)2 

4 (65)' 



81 (g2-r)(65 + r-g 



2^ 




rrDH,qq+, ^ ^ ( {Q' - Tf + {j^S + T - Q^f {Q^ - Tf + (65)^ 

(g2 - Tf 



81 65(65 + T - g2 



TL 27 (g2 _ T)' 

^ 2 (g2 - T)' + (65 + T - g2)2 

27 (^2^)' 



rrDH,qq'^, 2 (g2 - Tf + (65)' 
^TL ~ 



27 (65 + T-g2 



7.4 Soft-hard results 



(72) 



Here we hst the hard parts H{Ta) of the invariant projections T^^, T^^, T^^ and of the 
hadronic tensor, defined in section 2.2, for soft-hard processes. The results can be projected into 
a specific frame via eq. ( ^0|) and the matrices given in Appendix 7.1. This is mathematically 
trivial but it complicates the expressions substantially. For this reason we give only the invariant 



27 



projections. Note that especially i^TL = H{Tti)/2 for every frame described in section 2.2. 



""^'^^ ^~ 27 x,{Q'-T) 

™+,-. ^ ^ / ^2S + T-Q^ x^Q^S _ xg x,S{{x, + Xs)S + U-Q'') 
^ 27\{x, + Xs)S + U -Q^Q^-Uy xa-x, [Q^-Uf 

Xg {x, + x,){{x, - Xs){Q^ - T) + UQ' -U)- Xc^2S) - {x, - x,)Q^ 

Xa 6-5^ {Xc - - T) + 5(g2 _ ^) _ Q2 _ ^^^^^ 



27 \xa - X, (g2 _ T)(g2 _ T) 6-5 + T - g2 

x„s (xc + x,)(e2-5 + T - g2) + e2((xc + Xs)S + U-Q^)-Q^ + XcC2S 



XciQ"" -U) ^2S + T-Q 

H(m+') = -( ^ 



Xc^2S 2Xa{^2S + T-Q') ^ (Xc - Xa) (Q^ - T) 



HiTll- 



-9\ 



){Q^-T) ^2S{X,-Xa) X,^2S 

l Xa{C2S + T-Q^) 
4 .Te(g2 - T) 



JJ(^19+9^ ^ ( XgQ'^ {Xc + Xs)S + U 



4 \x,{Q^ - U) UQ^ - U) 
Xg (xe + Xs)S + U-Q^x,- Xs){Q^ - T) + 6(g' -U)-Q^- x,^2S 

xgx, {{xc + xs)s + u- g^)(6^ + T-Q^) + (g^ - T)(g^ - ^) - q^s 
xa-xc UQ^-uy 

ff(T^9+9^ ^ 1 / Xg ^,S + T-Q'{x, + Xs){Q'-T)-^2{Q'-U)-Q' + Xc^2S 
Xg (x, + Xs)(^2S + T-Q-')+ ^2{{xc + Xs)S + U-Q'') + Q^- x,^2S 

Xc UQ' - u) 

H(T'^3+^) = - ( - ^Q'iQ'~T)ixc-Xa)^2S + T-Q^ 

^ ^ 18 ^2-5 + T - g2 xli2S{^2S + T-Q^) ^28 

ff^'J'19+9'^ _ "^Xg* (Xn — Xn) 



9 X 



1 g2 f(xc + Xs)S + U-Q^ XsS s 



UQ'-U) {Q^-U)i2S + T-Q^ 

{XcS + U- g")((x, + Xs){^2S + T-Q^)+ ^2{{Xc + Xs)S + U-Q^)+Q^- Xc^2S) 

xMQ'-U)(^2S + T-Q^) 

2x„Q'((,r,, + .r:)S + U -Q') 
xlUQ'-U){^2S + T-Q^) 

TJ(rr19+9^ _ 1 ( '2xMQ'' -U)+ XcjXe + Xs){Q^ - T) + X%S + 2x2g^ , 

Ujxc + ^.j'g + ^ - g')(2g' + xciQ^ - T)) + 33:eg^(g^ - r) + 2g^ + xi^ls^ ' 

{Q'-U){i2S + T-Q^) 
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